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Abstract 



For quantum integrable models with elliptic i?-matrix, we construct the Baxter 
Q-operator in infinite-dimensional representations of the algebra of observables. 



1 Introduction 

The method of Q-operator, first proposed by Baxter in his solution of the 8-vertex model 
and the spin-| XYZ chain [|l|, is perhaps the most powerful approach to quantum in- 
tegrable models known today. As compared with the algebraic Bethe ansatz 0, one or 
another version of this approach, known also as the method of vacuum vectors, is in- 
dispensable for models which do not possess bare vacuum states. Among later striking 
developments we mention the works ||^, where the Baxter method was applied to the chiral 
Potts model. 

At that time, however, the Q-operator itself seemed to be a kind of virtual object 
for its explicit form was not known. In Baxter's original work, the Q-operator played an 
auxiliary role since the local Hamiltonian to be diagonalized was made of the quantum 
transfer matrix T(A). The main properties of the Q-operator Q(A) are the commutativity 
[Q(A), Q(/i)] = 0, [Q(A), T(/i)] = for all A, /i, and the following T-Q-relation: 

T(A)Q(A) = a(A)Q(A + 77) + 6(A)Q(A - r,) . (1.1) 

Here a(A), 6(A) are c-number functions depending on a particular model. Combined with 
(mo del- dependent) analytical properties of the Q(A), this relation allows one to obtain a 
system of Bethe equations and eventually to solve the diagonalization problem. In the 
absense of global bare vacuum states, Baxter's method to derive ( |1 . l]) was to consider a 
continuous-parametric family of special states that factorize into products of local vacuum 
vectors. Such states have especially simple properties under action of T(A). Treating 
them as columns or rows of some matrix Q(A), one arrives, after a proper normalization, 
at (|1 . 1|) . Some years ago this approach was extended to the case of twisted boundary 
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conditions. It remained obscure, however, whether the Q-matrix could play any other role 
than just a technical tool for deriving Bethe equations. 

The situation changed only quite recently, when it was realized that in integrable 
lattice models with an infinite-dimensional Hilbert space at each site it is Q rather than 
T that is the generating function of local Hamiltonians. With this motivation in mind, 
a considerable progress in constructing the Q-operator for a number of integrable models 
and in understanding its nature was achieved. 

The work |Q was perhaps the first, where the Baxter approach was adopted to inte- 
grable models with an infinite-dimensional quantum space. Specifically, in that paper the 
Q-operator for the quantum Toda chain was constructed as certain integral operator with 
explicitly known kernel. Recently, this remarkable development was extended to the 
integrable homogeneous XXX spin chain in a unitary discrete series representation of the 
non-compact group SU{1,1) (see also 0). 

For integrable models of XXZ-type, the nature of Q-operator was clarified in 0,0, 
where it was given a representation theory meaning. Specifically, the Q-operator can be 
obtained from the universal i?-matrix for the quantum affine algebra Uq{sl{2)), by taking 
trace of a chain product of the universal i?-matrices in a special g-oscillator representation 
in one of the spaces. The commutativity of the Q-operators, which comes out as a sur- 
prize in Baxter's approach, is now guaranteed by characteristic properties of the universal 
i?-matrix. The T-Q relation, from this point of view, is a simple consequence of fusion 
rules for representations of the quantum affine algebra. In fact, two different commuting 
Q-operators, Q'''(A) and Q^(A), can be introduced this way. They are two linearly inde- 
pendent operator solutions to the T-Q-relation ( |1.1| ) treated as a second order difference 
equation in A. As such, they obey a quantum Wronskian condition, which turns out to be 
even more instructive than the T-Q relation itself. 



As for matters more customary for theory of solitons, it was argued in |Ty, |Tl| that 
Q-operator is a quantum analog of Backlund transformation. In other words, it defines an 
integrable quantum model with discrete time by setting Q(A) to be the time shift operator. 
For any observable O, one sets 0(t + 6) = Q(A)0(t)(Q(A))^"'^, where the discrete time shift 
(5 is a function of A. Integrability implies that these discrete evolution operators commute 
for all 6. 

However, all this was about models associated with rational or trigonometric solutions 
of the Yang-Baxter equation (i?- matrices) . General i?-matrices with additive spectral 
parameter are parametrized via elliptic functions. The symplest elliptic i?-matrix is 

R(^) = — a rT~\ — ^a®c^a, (1.2) 



a=0 



where 6a{z\T) are Jacobi ^-functions, aa are Pauli matrices, and Uo is the unit matrix. 
This i?-matrix is associated with the 8- vertex model, being the matrix of local Boltzmann 
weights at a vertex. The transfer matrix of this model is the generating function of 
conserved quantities for the integrable XYZ spin-^ chain. By models of the XYZ-tjpe we 
understand more general integrable models with elliptic i?-matrix, for instance integrable 
anisotropic spin chains with higher spin. Integrable spin chains of XYZ-tjpe can be solved 
by the generalized algebraic Bethe ansatz [§, jT^ . 

However, the quantum affine algebra approach is not applicable in the XY Z-ca.se since 
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nothing similar to the universal elliptic i?-matrix is known. Nevertheless, the elliptic 
dependence on the spectral parameter does not destroy the functional relations between 
transfer matrices and Q-operators, including the T-Q relation and the Wronskian condition 
T3| , [0]. This suggests that everything what is possible to do in XXX and XX Z cases, 



might be possible, mutatis mutandis, for XYZ, too. 

In this paper we present a direct construction of the Q-operator for homogeneous A^- 
site spin chains of the XYZ-type. In lattice integrable models with elliptic i?-matrix, the 



algebra of local observables is the Sklyanin algebra |jT5|, [T^. A concrete model is defined 
by fixing a particular representation of this algebra. Similar to the s/(2)-case (models of 
XXX-type), the representations are labeled by a continuous parameter which is called 
spin, by analogy with s/(2). We allow the spin to take any complex value, so we work in 
an infinite-dimensional representation of the algebra of observables. In other words, we 
deal with representations of what would be called a non-compact real form of the Sklyanin 
algebra. Such representations can be realized by difference operators, so our Q-operator 
is a diffeference operator in many variables. 

Our strategy closely follows the original Baxter argument. First, we construct a family 
of local (infinite-dimensional) vacuum vectors that obey the celebrated "pair propagation 
through a vertex" property. Next, their product over all sites of the chain is, up to normal- 
ization, the kernel of the Q-operator. At last, one should prove the commutativity of the 
Q-operators. However, the practical realization of these ideas in the infinite-dimensional 
setting is by no means obvious. Technically, it is very different from what is customary 
in the 8-vertex model. The construction extensively uses such truly special functions as 
elliptic generalizations of the gamma-function and basic hypergeometric functions, and 
highly non-trivial identities for them. For example, commutativity of the Q-operators 
relies on an elliptic version of the Bailey transformation formula for the very-well-poised 
basic hypergeometric series loWg. The appearence of such a sophisticated mathematics in 
integrable models is not so surprizing. Recently, the elliptic series lo^Fg was shown |jl7| 



to offer a natural representation for "elliptic 6-j-symbols" , which are nothing else than 
solutions of the Yang-Baxter equation in a special basis. 

It is now a challenge to reveal a representation theory meaning of the elliptic Q- 
operator. Hopefully, this would make its properties and their proofs obvious. 

The paper is organized as follows. In the next section we present the main result in 
the shortest possible form. Sec. 3 contains the necessary material on rerpesentations of the 
Sklyanin algebra. The construction of the Q-operator begins in Sec. 4, where we find a 
family of local vacuum states for the quantum L-operator. In Sec. 5 we investigate direct 
products of such states along the chain and define unnormalized "pre-Q-operators" , Ql 
and Qr, which do not yet commute. Sec. 6 is rather technical. It is devoted to the proof 
of the basic commutation relation for the pre-Q-operators, QL(A)Q/j(/i) = QL(/i)Q_R(A), 
which entails commutativity of the Q-operators introduced in Sec. 7. In Sec. 8 we discuss 
more specific issues: a) twisted boundary conditions, b) the exceptional case of a "spin 
chain" with only one site, N = 1, and c) the trigonometric (XXZ-type) degeneration. A 
few concluding remarks are given in Sec. 9. 

Appendix A contains necessary information on the special functions involved in the 
construction. In Appendix B we give a direct proof of the intertwining relation between 
representations of spins i and —i — 1. 
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2 The main result 



To state the result, we recall the definition of the quantum transfer matrix T(A) for the 
XY Z spin chain with spin I. The elliptic quantum L-operator is the matrix 



L(A) = i 



/ ^i(2A)so + ^4(2A)s3 ^2(2A)si + ^3(2A)s2 
1^ ^2(2A)si - ^^3(2A)S2 ^i(2A)so - ^4(2A)s3 



(2.1) 



with non-commutative matrix elements. Specifically, are difference operators in a com- 
plex variable z: 

_ e,^,i2z-2iv) _ g,+i(-2^-2£77) 

9^{2z) 9,{2z) ' ^ ' 

introduced by Sklyanin |l6l. Here 9a{z) = 9a{z\T) are Jacobi 6'-functions with the elliptic 
module r, Imr > 0, £ is a complex number (the spin), and G C is a parameter which is 
assumed to belong to the fundamental parallelogram with vertices 0, 1, r, 1 + r, and to 
be incommensurate with 1, r. Definitions and transformation properties of the ^-functions 
are listed in Appendix A. 

12 N 

Let us consider the matrix product of L-operators (|2.1|) : L (A) L (A) . . . L (A), > 1, 

i 

where each operator matrix L (A) has the form (|2.1j ) with acting on its own variable Zi. 

i j 

Clearly, matrix elements of L (A) and L (A) commute for i j. The trace of this product 
(the sum of two diagonal elements), 

T(A)=tr[L(A)L(A) ... L (A)] , (2.3) 

is called the quantum transfer matrix of the spin chain. 

Our main result is the explicit realization of the Q-operator for this spin chain: 

N jv 

Q(A) = r-^(2A - 2ir]) ; H q.,^i(A): e^'+'^^^.-^^^^ , (2.4) 

i=l 

where 



^ T{2zi)T{zi-z,+, + 2X)T{z, + z,+, + 2\) ^ 

"^''^'^ ^ T{2zi + 2X + 2£7])T{z,-Zi+i-2er])T{z, + z,+i-2er]) ^ 

^7] 7] 2ri 2r] 



(2.5) 



?7,r 



(it is implied that -Zat+i = Zi). Here T{z) = r{z\T,2r]) is the elliptic gamma-function 
1181,1191: 



oo ^_g2^i((fe+i)r+2(m+i)„-^) / oo^ sin(7r A;(2;z - 2r7 - t) ) 

T{z\T,2r]) = 11 ,_o„i(kr+2mn+z) = exp 



k,m=0 



l_^2ni{kT+2mv+z) 2k sm{2TTkr]) sm{TrkT) ^ 
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and 6^/5(01; 0^4, . . . , ag; z |?7, r) is the elliptic analog of the very-well-poised hypergeometric 
series (with the operator argument): 



A:=0 



z 



where = Oj^o ^i(2«'7 + 2j77). For more information on the special functions F, gW^s 
see Appendix A. The normal ordering means that all operators e~^^^^j are moved to the 
right. 

We prove the commutativity [Q(A), T(/i)] = [Q(A), Q(/i)] = 0, and the T-Q relation 
T(A)Q(A) = (2A - 2£r/)Q(A + 77) + (2A + 2£r/)Q(A - i]) . (2.6) 

Another solution to this operator difference equation has a similar form (see ( |7.8| )). 

The limiting case r — 0, when elliptic functions degenerate to trigonometric ones, 
corresponds to the XXZ spin chain. The structure of the formulas ( |2.4| ), ( |2.5| ) remains 
the same while the special functions get simplified. The elliptic gamma-function becomes 
the g-gamma function (or the g-exponential function) while the elliptic hypergeometric 
series gW^s turns into the little g-Jacobi function represented by the basic hypergeometric 
series 20i- 



3 Representations of the Sklyanin algebra 

The aim of this section is to give the necessary preliminaries on representations of the 
Sklyanin algebra. We begin with a few formulas related to the quantum L-operator with 
elliptic dependence on the spectral parameter. 

Consider the quantum L-operator (|2.1| ), where are realized as in ( |2.2|) . The four 



operators obey the commutation relations of the Sklyanin algebraQ: 



(3.1) 



with the structure constants lab = ^^a+i(0)6'b+i(2?7). Here a,b = 0, . . . , 3 and {a,P,-y} 
stands for any cyclic permutation of {1,2,3}. The relations of the Sklyanin algebra are 
equivalent to the condition that the L-operator satisfies the "RLL = LLR" relation with 
the elliptic i?-matrix (|1.2| ). 

Plugging the difference operators ( |2.2| ) into ( |2.1| ), one can represent the L-operator in 



a factorized form, which is especially convenient for computations: 

L(A) =^i(2A+):y-^(A+,z) [ J^,^ ) V{X.,z):. (3.2) 



^The standard generators of the Sklyanin algebra [^5|, Sa, are related to ours as follows: Sa 
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Hereafter we use the convenient notation A-t = A ± fr], 6a{z) = 6a{z\^). The dots mean 
the normal ordering when the shift operators are moved to the right and 

Note that 

L(-A + i(l + r)) = _e— i-+4-'ALT(A) (3.3) 

where L""" means transposition in the two-dimensional space. 

The parameter i in ( p.2|) is called the spin of the representation. If necessary, we 
write Sa = s^^ to indicate the dependence on i. When i E these operators have a 

finite-dimensional invariant subspace, namely, the space B^^ of even 6'-functions of order 4i 
(see Appendix A). This is the representation space of the {2i + l)-dimensional irreducible 
representation (of series a)) of the Sklyanin algebra. For example, at £ = | the functions 
(^a{z), (^s{z) form a basis in Bg", and the generators s^, with respect to this basis, are 
represented by 2x2 matrices {—iY"'^ {9a+i{ri))~^aa- In this case, L(A) = R(A — ^rj), where 
R is the 8- vertex model i?-matrix ( |1.2|) . 

As is proved in [EH], the space B|^ for i G is annihilated by the operator 



2i+l 

w.= E(-i)' 

k=0 



2i + l 
k 



Uf=V0i{2z + 2{j-k)r]) 



lj=0 

Here and below we use the "elliptic factorial" and "elliptic binomial" notation: 

W! 



[j] = ^i(2jr/), M! = n[j] 



n 
m 



[m] ! [n — m] ! 



(3.5) 



Another important property of the operator established in is that intertwines 
representations of spin £ and of spin — (£ + 1): 



s 



= a = 0,...,3. (3.6) 



The same intertwining relation can be written for the quantum L-operator (p.l|) . The 
operator is an elliptic analog of [d/dzY^^^ in the following sense. In the case of the 
algebra s/2, the intertwining operator between representations of spins £ and —l—l (realized 
by differential operators in z) is {d/ dzY^^^ . It annihilates the linear space of polynomials 
of degree < 2£. 

The operator W£ can be extended to arbitrary complex values of £ in which case it is 
represented by a half-infinite series in the shift operator e^''^^. The series is an elliptic 
analog of the very- well-poised basic hypergeometric series with an operator argument. For 
its explicit form see ( [7.131 ). The intertwining relations ( |3.(i| ) hold true in this more general 
case, too. The proof is outlined in Appendix B. 

In what follows we need infinite dimensional representations of the Sklyanin algebra in 
some functional space. Operators ( p.2|) provide such a representation. However, they are 
not sufficient for our purpose since it is necessary to specify the functional space where 
these operators are going to act. It is not clear how to do that for continuous functions. 
On the other hand, the difference character of the operators ( |2.2|) allows one to restrict 
them to discontinuous functions of the form described below. 
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Let 6{z) be the function equal to zero everywhere but at z = 0, where it equals 1: 
6{z) = 0, z 0, 5(0) = 1. Clearly, z6{z) = 0. The Sklyanin algebra realized as in ( p.2|) 
acts, for any G C, in the space of functions of the form 

f{z) = Y.h6{z-iy + 2kr]), heC. (3.7) 

This space is isomorphic to the linear space of sequences {fk}kez- We call functions of the 
form ( |3.7|) combs. 

A comb is said to be finite from the right (respectively, from the left) if there exists 
M G Z such that fk = as k > M (respectively, k < M). Let (respectively, C^) be the 
space of combs finite from the left (respectively, from the right). 

We define the pairing 

{F{z), 6{z-a)) = F{a) (3.8) 
for any function F{z), not necessarily of the form (|3.7|) . In particular, 

{6{z-a), 6{z-b)) = 6{a-b). (3.9) 

Formally, this pairing can be written as an integral: 

{F{z), 6{z -a)) = J F{z)6{z - a)dz 

(perhaps a g-integral symbol would be more appropriate) . We stress that the integral here 
means nothing more than another notation for the pairing, especially convenient in case 
of many variables. By linearity, the pairing can be extended to the whole space of combs. 
We note that the pairing between the spaces and is well defined since the sum is 
always finite. 

Combs are to be viewed as kernels of difference operators. By a difference operator in 
one variable we mean any expression of the form 

D = ^ Cfc(z)e('^+2fcr?)a. ^ ^ ^ _ ^3^^q) 
fcez 

The comb 

D{z, = E ^k{z)5{z - C + /i + ^k'n) , (3.11) 
fcez 

regarded as a function of any one of the variables z, is the kernel of this difference 
operator in the following sense. Using the pairing introduced above, we can write: 

(D/)(z) = / D{z, OfiOdC = E ^k{z)f{z + /i + 2kv) . (3.12) 

The kernel D{z,Q can be viewed as an infinite matrix with continuously numbered rows 
{z) and columns {Q. Then the convolution with respect to the second argument of the 
kernel, as in ( |3.12|) , defines action of the operator from the left. The convolution with 
respect to the first argument defines the action from the right, 

(/D)(z)= J f{OD{C,z)dC, (3.13) 
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equivalent to the action of the transposed difference operator from the left: 

= ^ e-(M+2fc„)9.^^^^) _ ^ ^^^^ _ ^ _ 2A;77)e-('^+2'='')^^ . (3.14) 

The transposition t is the anti-automorphism of the algebra of difference operators such 
that In terms of the above pairing we can write (/, Dg) = (D*/, g). 



The kernels of Sklyanin's operators (^]2 



are: 



Let us also note that s^t = -{-lY-'^>ei{2z)si-^-^\9i{2z))-\ 




4 Local vacuum vectors 

We are going to adopt the Baxter method of vacuum vectors to the infinite-dimensional 
representations of the Sklyanin algebra described in the previous section. 

Let us recall the general definition of vacuum vectors. Consider an arbitrary L-operator 
L with two-dimensional auxiliary space C^, i.e., an arbitrary 2x2 operator- valued matrix 
whose matrix elements are operators Ln, L12, L21, L22: 

(4.1) 

They act in a linear space Ti which is called the quantum space of the L-operator. 

For a moment, let ip, denote vectors from C^. Let X G ?i, then acting by the 
quantum L-operator on the tensor product X ® 0, we, generally speaking, obtain a mixed 
state in the quantum space: LX®0 = Xi (8>0i -I-X2 ®02- The special case of a pure state, 

LX ® = X' ® , (4.2) 

is of particular importance. Given an L-operator, Baxter's receipe is to find all pairs X, 
such that ( [4.2| ) holds. Taking the scalar product with the vector ijj-^ orthogonal to ip, we 
get: 

(^^L0)X = O, (4.3) 

i.e., the operator K = ('0"'"L0) (acting in the quantum space only) has a zero mode X ETi. 
Suppose ( [4.2|) (or ( [4.3|) ) holds with some vectors 0, -0; then the vector X is called a vacuum 
vector of the L-operator. 

A few remarks are in order. The relation ( [4.2|) (in the particular case = C^) was the 
starting point for Baxter in his solution of the 8- vertex and XYZ models IQ. This is what 
he called the "pair-propagation through a vertex" property. In the context of the quantum 
inverse scattering method p| the equivalent condition (|4.3|) is more customary. It defines 
local vacua of the (gauge-transformed) L-operator. A generalization of that solution to the 



higher spin XYZ model was given by Takebe in [|1^. An algebro-geometric approach to 



the equation ( |4.2| ) for finite-dimensional matrices Lj^ was developed by Krichever 
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For L-operators with elliptic spectral parameter it is convenient to pass to the elliptic 
parametrization of the components of the vectors 0, tp. We introduce the 2-component 
vector 

and the corresponding covector {(\ = (^4(C)) ^3(0) (recall that ^a(C) = ^a(Cli))- The 

^'^(C) )' '^^^^'^ ~ ^' generally, we have: 

(e|C)^ = 201(^ + 0^1(^-0- (4.4) 

Note also that 

\C + hil + r)) = e-^-'-^^\0^. (4.5) 

Let us specify the general definition of the vacuum vectors for the elliptic L-operator 
( p.l|) . Writing L(A)|0 (respectively, (OL(A)) we mean that the 2x2 matrix ( p.l|) acts on 
the 2-component vector from the left (respectively, on the 2-component covector from the 
right), as usual. Similarly, we introduce right and left vacuum vectors Xr, Xl according 
to the relations 

(CIL(A)X^ = {^\X'^ , X^(C|L(A) = Xi(e| . (4.6) 
In the latter formula the difference operators act on X^ from the right. 
The following pictorial representation is helpful: 



(CI 



X 



■10- 



(0 



Here we use the standard conventions of the quantum inverse scattering method. 
Introducing the operator 

K = K(C,0 = (C|L(A)|0\ (4.7) 



we can rewrite ([4.6|) as KX/j = X^K = 0. The explicit form of the operator K is easily 
found from (|3^): 

K = K(C, = P(Oe^'^ + p{-z)e-'^'^ , (4.8) 

where 

P(0 = Ujii^ + < - KWliz + + A_) . 

So, the equations for the right and left vacuum vectors read 

p{z)Xr{z + r]) = -p{-z)Xn{z - r]) , (4.9) 

p{-z - + v) = -p{z - v)Xl{z - T]) . (4.10) 

These equations themselves are not yet enough. An important additional input is the 
space where to solve them. Our choice, justified by the final results, is the following. We 
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are going to solve the first equation in the space of combs finite from the left, , while 
the second equation is to be solved in the space of combs finite from the right, C^. (For 
precise definitions of these spaces see Sec. 3.) 

Let us begin by solving equations (|4.9|) , ( |4.10|) in the space of meromorphic functions 



of z. This is easy to do in terms of the elliptic gamma-function r(2:|r, 2r]) introduced by 
Ruijsenaars [0 and recently studied in detail by Felder and Varchenko |]19[. This function 



was introduced just to solve equations of such a type. Its main properties are summarized 
in Appendix A. In what follows we write r{z\T,2r]) = T{z) for brevity. Let Xii{z), Xl{z) 
be meromorphic solutions of eqs. ( |4.9|) , ([4.10|) , respectively. Using (|A9|) , we easily find: 

T[z + eC - A+ + r])V{z + + A_ + ?]) 
X^{z) = e'-'^'+'^'M^\)9,{2z) n ^^'^f'^-'ll^'^'i^^-^ r . (4.12) 



The 2;-independent normalization factors C(j(A,C — i'^ ( [4.11| ) and $£(2A) in ( [4.12| ) are 
introduced here for later convenience. The function c<j(A, Q is given by 

, eii2\-2M2ri) I niC „A . ^ 

(the ^-functions here have the modular parameter 21] rather than r!). The function $^ 
obeys the difference equation 9i{z - 2£ri)^i{z + 2r]) = 6i{z + 2{i+l)r])^i{z). For £ e 
the solution is very simple: 

£ 

^ei^) = n di{z + 2kri), (4.14) 

k=-e 

where k ranges over —i, For arbitrary complex values of i the function 

can be defined through the elliptic gamma-function (see ( |A16|) ). If necessary, we write 
XrA^) = ^rA^i C, ^, A). Note that XrAz; C, ^, A) differs from XrA^i C, -A) only by 
a ^-independent factor. 

An important remark is in order. The solutions ( [4.11| ), (|4.12| ) are by no means unique. 



Clearly, one can multiply them by an arbitrary function f{z) such that f{z + ri) = f{z — r]). 
In this paper, we do not take care of such 2?7-periodic factors fixing them to be constant. In 
other words, ( ^.11| ), ( [4.12[ ) are "minimal" solutions to eqs. ( |4.9| ), ( |4.1CI| ) in the sense of []TB|. 



By this we break the symmetry between the three periods 1, 2?7, r involved in the problem. 
In particular, the modular symmetry is broken since 277-periodic factors may be important 
for good modular properties of the solutions. Indeed, the modular transformation of the 
elliptic gamma- function has the form |T3: 



r(z|r,2r/) = e^"-^(") " 



p/ Z — T I _1_ T_\ 

\ 2rj \ 2r)'> 2r]' 



where P{z) is a polynomial of third degree (see ( A12 )). The denominator is a 2?7-periodic 



function which does not enter the minimal solution of the modular transformed equation. 



We note in passing that this formula is an elliptic version of eq. (40) from the paper p2 
where a modular double of quantum group was introduced. We hope to discuss this point 
elsewhere. See also the end of Sec. 8. 
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Now we can proceed to built a half-infinite comb solution. Any infinite comb of the 
form X{z) = Xji^l{z) J2kez^i^ — zq + 2kr]) automatically satisfies eq. (resp., (K:.10|) ). 
How to truncate the comb holding these equations valid? Let X{z) be any one of Xr{z), 
Xl{z). The position of zeros and poles of the F-function suggests two possibilities: 

(1) To set Zq to be a zero of the function X: X{zq) = 0. If this zero comes from a pole 
of one of the F-functions downstairs, then X{zq — 2kri) = 0, k > 0, and the comb 
is truncating from the left. If zq is a zero of one of the F-functions upstairs, then 
X{zo — 2kri) = 0, k < 0, and the comb is truncating from the right. 

(2) To set zq to be a pole of X and define X{z) = J2k(^z X{z)6{z — Zq + 2kri). 
If Zq is a pole of a F-function upstairs, then the comb is truncating from the right. 
If the pole comes from a zero of one of the F-functions downstairs, then the comb is 
truncating from the left. 

In the sequel we use both possibilities: (1) for Xr and (2) for Xl. Other options are to 
be further analysed. 

Assume that the parameters are in general position, i.e., there is no cancellation of zeros 
and poles in ( [4.11|) . (For £ G |Z+ the sufficient conditions are A 7^ 2m?7, ^ — C 7^ 2m?7, 



m G Z.) Then — — rj — 2kr]) = for all integer k > 0. At the same time this 

function is not identically zero when k is negative. Therefore, 

00 

X+(z; C, A) = Xr{z; (, ^, A) J] 5(z - ^ + A_ + r/ - 2kT]) (4.15) 

k=l 

is a non-trivial comb finite from the left which solves eq. ([4.9| ) . Clearly, this is not the only 
possibility to truncate the infinite comb from the left, even if to stick to (1). There are 
three other series of zeros of the function ( ^4.11| ) coming from poles of the gamma-functions 



in the denominator. Any one of them can be used to generate, in the same way, a half- 
infinite comb solution that belongs to . However, inasmuch as the vacuum vectors are 
going to be kernels of difference operators, we want the 5-functions in the vacuum vectors 
to depend on the differences z — ( or z — C,- So, along with (|4.15| ), we have another comb 
obeying eq. ([4.9|): 

00 

Xiiz; C, e, A) = Xniz; ^, C, -A) ^ 5(;^ - C - A+ + r/ - 2fcr/) . (4.16) 

k=i 

For Xl we use (2): Xl{z) has poles at z = C + A+ — 2kri, k > 0. In general position, 
these poles are all simple. The lattice of poles is finite from the right. Therefore, we can 
construct a solution to eq. ( [4.10|) in the space of combs finite from the right as follows: 



X+(A,C,e;^) = Eres M^'X,^, X) S{z - ( - X+ + 2kv) (4.17) 

k=0 ^-C+A+-2fc»? 

Similarly to the case of the right vacuum vectors, we have a second solution: 

00 

X^(A,C,e;^) = Eres| A'4^;e,C,-A)5(^-e + A- + 2fcr/). (4.18) 

k=0 ^ ^ 
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To make these formulas closer to ( |4.15| ), ( |4.16| ), we introduce the auxiliary function 



r(^ - C - A+ + e)T{z + C - A+)r(^ - ^ + A_ + e)T{z + ^ + A_) 
^ r(z - C + A+ + 2ri)T{z + C + A+ + 2ri)T{z - ^ - A_ + 2ri)T{z + ^ - A_ + 2r/) ' 
then, up to an irrelevant constant factor, 

X+iz; C, A) = lim eXl{z; (, ^, A) J] - ( - A+ + 2k7^) , (4.19) 

k=0 
oo 

(A, C, z) = lim eXliz- ^, (, -A) ^ - ^ + A_ + 2kr]) . (4.20) 

fc=0 

Here one should first substitute the values of z in the 5-functions and then take the limit. 

Note that the function p{z) and the meromorphic solutions A/j^ are invariant with 
respect to the simultaneous change A —A, C ^ ^ while the half-infinite combs are 
transformed one into another: 

X^(2;C,e,A) = X+(^;e,C,-A), 

(4.21) 

The construction of local vacuum vectors is completed. The vacuum vectors are de- 
fined up to an arbitrary normalization factor which may depend on and A. In our 
normalization, the final formulas have the most simple form. 

By a straightforward computation, which uses some standard identities for theta- 
functions, one can prove the following relations valid for = X^^l: 

{CU\)Xr{z- C, A) = e,{2\_mXR{z- C, A + r^) , (4.22) 

Xi(C, A; z)(C|L(A) = 0i(2A+)Xi(C, iA-W. , (4.23) 



which are simply ( |4.6| ) in the specified parametrization. In the first equation, each operator 
entry of the matrix L(A) acts on the vacuum vector from the left while in the second one 
they act from the right. In the both equations, the two-component covector {(\ is applied 
from the left. Using ( |3.3| ) and ( [4.5| ), one can represent these equations in the form, where 
the two-component vector is applied from the right: 

L{X)\O^Xn{z; C, A) = ^^i(2A+)|C)^X^(^; C, ^, A - r/) , (4.24) 

X^(C, A; ;^)L(A)|0^ = ^i(2A_)|C)^Xz.(C, ^A + V;z). (4.25) 



For the holomorphic vacuum vectors ( (|4.11|) at £ G Z_|_) these relations hold [12] under the 
additional restriction ^ — C G 2r]Z. 
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5 Global vacuum vectors 



The next step, in accordance with [Q], is to find the vacuum vectors for each L (A) and 
consider the action of the transfer matrix T(A) (|2.3| ) on their direct product along the 
chain. It then appears that the result of this action has a very special and simple form: 
it is a sum of two terms, each of them being a vacuum vector that belongs to the same 
family. From the technical point of view, the easiest way to see this is to use the gauge 
transformation argument, now standard. 

Let us make the "gauge transformation" of the L-operator: 



L(A) ^G(0)L(A)(G(Cm))-'=L' (A) 



(5.1) 



where G{() are some 2x2 c-number matrices depending on the parameter (. If Ca^+i = Ci, 
then the trace is not changed: 

/I 2 N \ 

T(A)=tr L' (A) L' (A)... L' (A) . 



We set the matrix G{() to be 



^(0 



1 



(5.2) 



Then, acting on a vacuum vector Xji by each element of the gauge-transformed L-operator 
and using the relations ( [4. 221) , ( [4.24|) , we obtain: 



^2iXr L122XR 



1 




-^"3(0 
^4(0 







^"4(0 ^„ 



(5.3) 



Here X'j^ = ei{2\_)XR{zi; C,^, A + r]), X'^ = ^i(2A+)Xfi(z,; C^, A - v) and we imply that 
Q = (, Q^i = ^. The star in the left-down corner stands for something irrelevant for us. 

We define the global vacuum vectors X^j^ to be products of the local ones over all sites 
of the chain: 



N 



x±(A) = x^{z, c, A) = n 0, 0+1, A) , 

i=l 
N 

X±(A) = X±(A, C, z) = Y[ X±(A, Q, 0+1; , 



(5.4) 



(5.5) 



i=l 



where the notation z = {zi, . . . , z^}, ( = {Ci, • • • , Ca^} is used. The parameters Q are 
naturally asociated with the horizontal edges of the lattice: 
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Ci 


C2 


Cs 


Cn 


Ci 











Cn+1 — Cl 



zn 



Then 

T(A)X|(A) = tr n (l' (A)X|(z,; 0, O+i, A)) 

and similarly for the X^(A). The triangularity of the matrix ( p. 3D allows one to compute 
the trace. Taking into account the periodic boundary condition, we arrive at the following 
relations: 

T(A)X±(A) = (2A_)X|(A + r^) + <(2A+)X±(A - r/) , (5.6) 
Xt{\)T{\) = en2X-)XtiX + + <(2A+)X±(A - v) ■ (5.7) 

Let us introduce left and right Q-operators, and (the pre-Q- operators). They 
are difference operators acting on functions f{zi, . . . , zn) of variables, respectively from 
the left and from the right, as follows: 

(Ql/)(C) = / dzXiiX, C, z)f{z) , (5.8) 

(/Ql)(C) = / dzf{z}X^{z,CX). (5.9) 

The kernels here are the global vacuum vectors defined above. The symbol J dz = 
J ■ ■ ■ J dzi . . . dzN means the natural pairing in the tensor product of the local functional 
spaces. On each factor, this pairing coincides with ( |3.8| ). 

The formulas ( ^.6|) , ( ^.7|) can be written as the operator relations: 

T(A)Q±(A) = <(2A_)Q±(A + + <(2A+)Q±(A - v) , (5.10) 

Q±(A)T(A) = <(2A_)Q±(A + + <(2A+)Q±(A - v) ■ (5.11) 

One might multiply the pre-Q-operators by a scalar normalization factor /(A) such that 
/(A) = /(A + rj). Clearly, these relations remain unchanged. 

6 Basic commutation relations 

In this section we prove the basic commutation relations 

Ql(A)Q|(A') = Q^(A')Q|(A), (6.1) 

Q^(A)Q^(A') = Q^(A')Q^j(A). (6.2) 
Let us start from the proof of ( |6.1| ) for Q^^. To this end, consider the function 

J d'^zXtiX, C, z)XUz, C, A') = n y dz^XUX, 0, 0+1; z.)X^{zf, C, C'+i, A') (6.3) 



i=l 
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As is clear from ( |6.3| ), the kernel is a product of local multipliers associated with each pair 
+ of adjacent sites of the chain. For a while, we set Q = (, Ci+i = ^- Therefore, we 
should investigate the function 



(6.4) 



Let us begin with the X^. Up to an irrelevant constant factor, (|4.19|) reads 
Xt{\ C, z) = limfeA-KC + A+; C, A)] ^fc(C, A)5(^ - C - A+ + 2kr^) , (6.5) 



fc>0 



where 



X 



1 ^i(2C + 2A+-4^ 
xu[Q. A) - — ^^(2C + 2A+) 

'=-iei(-2A+ + 2jr/)^^i(-2A+-2C + 2jr/)^^i(^-C-2£r/ + 2jr7)ei(-e-C-2£r/+2jr/) 



n 



^i(-2C + 2(j + l)r/)^i(e-C-2A + 2(j + l)r/)^^i(-e-C-2A + 2(j + l)r7) 



(6.6) 

Substituting (|;D and X^(;z; C', A') = A'^j(z; A') Efc'>o ^(^-'^' + A'_ -r?-2A;'r7) (see 
( [4.15| )) into (|6l^) , performing the 2;- integral according to the pairing rule (|3.9|) , and passing 
to summation over n = + fc' > 0, we get: 



Fg,(A,A') = liml£A'|(C + A+;C,e,A) 



n 

X E 5(C-r+ A + A'-r/-2nr^) ^ ^^.(C, \)Xn{C + K-2kv; C, e, A') • 

n>0 A:=0 



(6.7) 



The next step is to identify the finite sum over k with the terminating modular hyper- 
geometric function 10 Wg. For the definition see Appendix A. Plugging here the explicit 
formula for X^f ( [4. Ill ) and using ( [A13| ), one can represent ( |6.7| ) in the form 



Fg(A,A') = 4|i^^'(A, AO u{X', C--a F<HX)Ff^'{X, X') : x 

X YioWg{ai; 04, ... , ttg, -n)5{C-^' + \ + \'-r]-2nr]) , 

n>0 



(6.8) 
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where 



sym Y{C-^ + 2{£+l)7])T{CH + 2{i+l)v) H r(C + + A + A' + r^) 
is a symmetric function of A, A' for all (, ^, uj{X', C~0 is defined in ([4.13|) , and 

<7CC'(A) = r(c + C' + A + r^)r(c -C + x + v), (6.9) 

^ ^ _ ,-4...A r(c-e + 2A)r(c + e + 2A) 
- " r(2A_)r(2C + 2A+) — ' ^^-^^^ 

pcc'r.^ _ r(C + r + A+ - A-_ + mc - ^ + A^ - A-_ + y) 

r(e + e + A-A' + ,)r(e-e + A-A' + ,) • 

The function 10W9 is defined in (|A20|) . It is an elliptic analog of the very- well poised 
hypergeometric series. The values of parameters are as follows: 

-2r]ai = 2C + 2X + 2er] -2r]a7 = C+^ + 2ir] 

-2r]a, = 2X + 2iri -2r^as = C-^+2iv 

-2r7«5 = C + C'+A-A'-r/ -2r]ag = C+^' +X + X' -r] ^""'^^^ 

-2r]aQ = C-C + X-X'-T] -2r]aio = C-^'+X + X'-r] = 2n7] 

Note that the parameters in the right column are symmetric under the permutation A A'. 
The series is balanced that means that the balancing condition 

7 

2 + 3ai = J2 (6.13) 

m=l 

is satisfied. The series is terminating since ctio = —n, where n is a non-negative integer 
number. 

Let us move the factor F2(A, A') inside the sum J2n>o &^d set ^' = C + A + A' — 77 — 2nri 
in the n-th term (this is what the 5-function suggests to do). After using ( |A13|) , (|A14|) 
and extracting the common factor 

F«rA A'^ - r(2C + 2A+)r(-2A; + 2ry) 

^ ^ - r(C + e + 2A)r(-C + ^-2A' + 2ry) ' ^^"^^^ 

the sum over n becomes equal to 

T?at\ x/N ["l-"7 + l]n[/5l-tt7-«8 + l]n ... / x./^ c' , \ , \' O ^ 

^3 (^) X ) 2^ ? -TT-j^ -7-1 ioWg{ai; ^4, . . . , -n)6{C-^ +A+A -ri-2nr]) , 

(6.15) 

where we use the elliptic Pochhammer symbol = + 1] . . . [x + n — 1], and (3i = 
2q;i + 1— 04— as— ag. Now, we recall the elliptic analog ([A21|) of the Bailey transformation 
formula for the terminating balanced loWg. It is easily seen from ( |6.12| ) that the Bailey 
transformation of the parameters —>■ jSi (|A22| ) is equivalent to the permutation A ^ A' 
while all the other parameters remain unchanged. Therefore, as it follows from ( |A2 1| ) , 
each term inside the sum over n in ( |6.15|) is symmetric under the permutation X ^ A'. 
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The common factor is not symmetric. Combining it with the remaining factors and 
taking into account ( |A10 ), we obtain: 



(6.16) 

X E-^lym '(A, A'; n) cu(A, (-0 ^(A', C -e')'5(C-e'+A + A'-r/-2nr7) , 

n>0 

where J^sym^'(A, A'; n) is symmetric under the permutation of A and A' for all n, (, ^, 
and the functions g, u are given in (|6.9| ), ( [4.13|) respectively. 



To obtain the kernel (|6.3|) , we should set C = 0; ^ = d+i and take the product of 
the local factors ( |6.16| ) over i from 1 to provided (n+i = Ci- Therefore, the ratios 
of the (^-functions cancel. The only factor which is not manifestly symmetric under the 
permutation A ^ A' comes from the product of the ti;-f unctions. This factor is equal to 



N \ N 



exp ( - 1^ E(C; - Ci+i)' ) n - 0+1 - 2A|27^) 9,{C - CUi - 2A'|2r/))"' . (6.17) 



i=l 



However, we should take into account the constraints imposed by the 5-functions: (■ 
Ci+i — A — A' + ?7 + 2nj_i?7, where rii G Z. In other words, ( |6.17D is zero unless (I — Cj'+i 



— Ci + 2(nj_i — rii)?]. From the monodromy properties of the 6'-function it then follows 
that this function is symmetric, too. Therefore, the commutation relation (|6.1|) for Q^^ 
is proved. To prove it for Q^j^, there is no need to repeat the calculations since X£ is 
related to X^j^ via (|]2l]). 

The proof of ( |6.2| ) requires a similar but different calculation. Its scheme is outlined 
below. The kernel of the operator in the left hand side is 



Y[ J dziX^X, Ci, Ci+i; Zi)Xj^ {zi] Ci, C'+i' A') — 

(6.18) 

n / dZiX^{\,(i,Ci+l] Zi)^RiZi'Xi+l,Ci, —^') , 



where ( |4.21| ) is used. Again, set Q = (, Ci+i = C,, Ci = C'+i = C and examine the 



symmetry properties of the function 

F^?|,(A,AO = J dzXt{XX,i;z)X^{z-iX,->^) 

under the simultaneous change C ^ ^) C ^ A — * —A', A' — > —A. It is enough to show 
that the product of such functions over the chain is invariant under this substitution in 
each factor. Proceeding as above, we obtain the representation of this function similar to 

I: 

Fg(A,AO = FlIi^^A, AO ..(-A-, ^-CO Ff^(A)FF'(A, \') : x 

(6.19) 

X Eio^9(«i; "4, • • • , aio)5(C-C'+A-A'-r7-2nr/) 

n>0 
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with the same functions -Fsym; Fi, F2, and g as before. The parameters of the series 
10 Wg are also the same (see ( |6.12 )); however, for the sake of convenience they are permuted 
in the following way: ai = ai, = ag, 0:5 = ag, = aio, ay = a^, ag = 05, ctg = a^, 
<^io = c^e- (Clearly, the series is symmetric under the permutation of all parameters except 
«!.) What differs is the argument of the ^-function. Note that the series is terminating 
since the 5-function is zero unless aio = —n. Next, we make the Bailey transformation 
( |A21|) passing from to Pi according to ( |A22| ). The result can be represented as a product 
of three factors: 



Fg(A,A')= 4|,(A,A') X 



invariant under 



factors that cancel in the 
product over the chain 



where 



44 (A, A') 



exp(-7r^(r-C0V(2r7)) 
9i{e-C-2X'\2vWi{C-^+2X\2r]) 



The product Hili A^]^'^^ (A, A') is the only thing here which is not manifestly invariant 



under the simultaneous change Q ^ Q+i, (■ ^ Ci+i? A 



-A' in each factor. However, 



the same argument as the one used for (|6.17|) establishes the invariance. 



7 From pre-Q-operators to Q-operator 

Suppose the operators and are invertible at some particular A = A^,. Following 
we define 

Q=^(A) = (Q^(A.))-^Q^(A) = Q|(A)(Q+(A.))-i . (7.1) 

The second equality is due to the basic commutation relations ( p.l|) , ( |6.2|) . It then 
follows, on applying (|0|), (p^ ) once again, that [Q+(A), Q'^{fi)] = [Q+(A), Q~(/i)] = 
[Q-(A), Q-(^)]=Ofor all A, /i: 

Q+(A)Q+(/i) = 

= (Qt(A.))-iQ+(A)Q+(/x)(Q+(A.))-i = 

(7.2) 

= Q+(/x)Q+(A), 



and similarly for Q . Both equations (|5.1CI|) , (|5.11|) lead to the Baxter T-Q relation 



T(A)Q(A) = (2A - 2£7])Q{X + r]) + (2A + 2£r/)Q(A - r]) , (7.3) 

where Q is any one of Q^. Besides, it follows that Q^(A) commute with T(/i) for all A, /i. 

Suppose i ^ Z+. Then the convenient choice of the normalization point A,, is A* = —£77. 
Indeed, at A = —ir] the series (|6.5|) terminates at the k = term; therefore, Q+(— £77) is a 
diagonal operator. Since we did not take care of ^-dependent constant factors, we continue 
to do so and redefine the Q-operator dividing it by the constant {—Air]). Note that this 
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normalization allows us to define the Q-operator for £ G Z+ as well. We express the result 
through the function (|A16 ): 



QU-^V) = n - + • (7.4) 

i=i 

Thus, by virtue of ( |5.8|) , the kernels of the Q-operators are 

TV 

Q±(A; C, ^) = n ?^(^; 0, Cm, , (7-5) 

i=l 

where 

+(. . . . r(20)r(o - 0+1 + 2A)r(o + c+i + 2A) 

^ r(2A_)r(20+2A+)r(c.-Cm-2^r^)r(0+Cm-2^^) 

(7.6) 

X XI ^fc(Ci'Ci+i, -^)^(2:i - - + 2fc?7) 

fc>0 

(xfc is defined in (|0|)), and g~(A, Q, Q+i, Zi) = g"'"(— A, Ci+i, 0? -^i)- Representing this in the 
operator form, according to the rule described in Sec. 3, we obtain: 

^ ^ ' T^i2X-2£r])l{T{2z, + 2X + 2£r])T{z,-Zi+,-2i7])T{zi + Zi+i-2ir]) 



X :U^W,{-'-^-i; -^-i, -'-i±^-i, e-^%|^,r): X 



X 



TV 

exp((A + £r/)X9 



Zm y 5 

m=l 



(7.7) 

1 r(2z,)r(^,-^,_i-2A)r(z, + z,_i-2A) 

^ ^ ' T^{-2X-2£r])l}^T{2z,-2X + 2er])T{z,-z,-i-2er])T{zi + z,-i-2£7]) 



N 

X 



:U^W,{-'-^-i; ^-i, -5l±^-£, e-^%|,,r): X 

jJl ^ 1] 1] 2r) 2rj ^ 



N 

X exp{i-X + iv)T.d.^)C, 

m=l 

(7.8) 

where C is the cyclic shift operator, 

(C/)(^l, Z2, . . . , Zn-1, Zn) = f{z2, Z3, . . . , Zn, Zi) , (7.9) 

with the kernel C{C,z) = Ylf=i5{zi — Ci+i)- Clearly, C commutes with the T(A) and the 
Q^(A). 
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In general, the series qW^ in (|777|), ( fTsP is not terminating. However, at A = ^iri±mri, 
m G the series qW^ in (|7.7| ) (respectively, in ( |7.8|) ) terminates at the m-th term. Note 
also that 

q^{±er]Tmri) = 0, m = 0,l,2,..., (7.10) 
due to zeros of the scalar common factor. At the points ±£1] we have: 

where I is the identity operator. So, at £ G the Q-operators have poles at these points. 

Let us evaluate the Q-operator at another distinguished point. Note that if A = (£+l)?7, 
then (i in the kernel of Q"*" decouples from Q+i, so the operator simplifies: 

N N . 

Q+((£ + l)r]) = (n M^i - Zi+i)Mzi + Zi+i)) n ■ (7-12) 
i=l j=l 

The operators W£ pairwise commute since they act on variables attached to different sites 
of the chain. From ( [7.7|) one obtains the representation of the operator through the 
elliptic hypergeometric series 4W3 with an operator argument: 

T{2z + 2{2i+l)r])' r] ' ' ^ * ^ ^ 

At £ G iZ+ the series is terminating and gives the operator introduced in ( p.4|) . So, at 
A = (£ + l)rj we obtain a distinguished factorized integral of motion 

N N . 

A = (n ^eizi + Zi+^)<^eiz, - Zi+^)) H w, (7.14) 
i=i j=i 

(no normal ordering is required since the we commute). 

Let us give some remarks on the operator A. It is clear from Sec. 3 that for i G 
|Z+ the operator A annihilates the space (O^)*^^. At £ = 0, this integral of mo- 
tion is proportional to T(0) and has especially simple form: 2^ (jl^i6i{zi + Zi+i)6i{zi — 

+i)/6i{2zi)^ Ylf=i sinh{ridz)- (In the limit to the XXX spin chain it coincides with the 
integral of motion found by Lipatov in the context of high energy QCD, see also ||2^ .) 
Presumably, for i G A can be represented as the trace of the product of more general 



z. 



L-operators obtained from (|2.1| ) by the fusion procedure, taken at a particular value of A. 
For brevity, we denote ipi{zi,Zi+i) = ^i{zi — Zi+i)^e{zi + Zi^i). Together with (|3^), the 
commutativity of A and T(A) implies the intertwining relation 



N N 



tW(A) n^K^.,^m) = (7.15) 



i=l i=l 



where the superscript of T indicates the value of spin. A direct proof of this relation is 
given in Appendix B. We stress that the intertwining relation ( p.6|) holds for all the spin 
operators separately while ( [7.15| ) is valid for the transfer matrix only. 
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The operators Q^(A) are two operator solutions of the Baxter T-Q relation ( [7. 3D treated 
as a difference equation in A. Therefore, their Wronskian, 

W(A) = Q+(A + r^)Q-(A) - Q+(A)Q-(A + r^) , 

obeys the functional equation (2A - 2£r])W{X) = 0f (2A + 2£r/)W(A - r]). Taking into 
account ( [7.10|) , one can write 

W(A) ^<(2A+2£,|2„) (£||±|i±|))"p(,)AC. (7.16) 

where p(A) is a diagonal operator such that p(A + 1) = p(A), p(A + rj) = e^'^*^^p(A), and 
p(£?7) = I. The operator A is given in ( |7.14| ). (The simplest way to find the operator 
factors is to calculate the Wronskian at A = irj.) To find p(A) explicitly, a more careful 
analysis is required. 



8 Miscellanea 

Twisted boundary conditions 

The twisted boundary condition means inserting a c-number matrix inside the trace of 
the chain of L-operators (|2.3|) . To preserve integrability, this matrix must be a constant 
solution to the RLL = LLR relation. In the elliptic case, there are only three such solutions: 
the Pauli matrices cTq, a = 1, 2, 3. Therefore, the definition of the quantum transfer matrix 
with twisted boundary conditions is as follows: 

T(A)=tr[L(A)L(A)... L (A)a«] . (8.1) 



For definitness, we consider the case of cti 



1 

1 



Again, we employ the gauge transformation ( p.lD with the same matrix ( p.2| ). The 
whole argument remains the same with the only modification that Ca^+i Ci- The Ca^+i 
should be adjusted in such a way that the matrix (t[ = G{(iy^i)(Ji{G{(i))^^ be lower- 
triangular. Calculating it explicitly, we obtain 



/ HCn+i) -202(Cl + CiV+l)^2(Cl-CjV+l) 

V 1 -UCi) 



(8.2) 



where a simple identity for 6'-functions has been used. Therefore, if Ca^+i = Ci + |) then 
the right upper element is zero. 

With this modification, the definitions of the vacuum vectors ( ^.4|) , (|5.5|) and of the pre- 
Q-operators (|5.8|), ( p.9| ) remain the same. However, equations ( |5.1C1| ), ( |5.11| ) get modified: 

T(A)Q±(A) = (2A_)Q|(A + r^)- <(2A+)Q±(A - v) , (8.3) 

Q±(A)T(A) = <(2A_)Q±(A + <(2A+)Q±(A - v) • (8.4) 
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The new constraint on Civ+i and Ci does not destroy the basic commutation relation 
( |6.1|) . Indeed, the (/-factors in (|6.16|) still cancel in the product over the chain since 

^CiC( = ^Civ+iCWi (El) ^"^^ 

The normalization of the Q-operators goes through as before. The T-Q relation is 

T(A)Q(A) = (2A - 2£ri)Q{\ + r]) - (2A + 2^r7)Q(A - r/) . (8.5) 

The Q-operators are given by the same formulas ( |7.7| ) ( [7.8| ), where it is implied that 
Cat+i = Ci + |) and the cyclic shift operator C should be replaced by the operator C such 
that {C'f){zi, 2:2, ... , zn) = f{z2, Z3,..., zn, zi + \). 

The case = 1 



For = 1, there is only one L-operator in the chain ( |2.3| ). Trace of the L(A) is the simplest 
quantum transfer matrix T(A). It is proportional to Sq: 

T(A) = trL(A) = 0,(2A)s„, (8.6) 



^1(22 + 2£r/) 



7) 



^1(22) eii2zj 

This is the difference analog of the Lame operator proposed in ||2^ . The difference operator 
given by eq. ( |8.7| ) for positive integer values of £ and arbitrary generic rj has 2^+1 stable 
bands (and 2£ gaps) in the spectrum. 

The construction of the previous sections gives a remarkable family of commuting 
difference operators in one variable. For £ G Z+ this family of operators commuting with 
( ^.71) coincides with the one found in |2^ by means of a different argument. 

The case = 1 is exceptional and must be considered separately. Some of the argu- 
ments for N > 2 are not directly applicable to = 1 since the vacuum vectors ( |4.15|) , 
( [4.17| ) are to be taken at ( = ^ that for £ G Z+ means the non- generic position of the pa- 
rameters. Nevertheless, it is a matter of direct verification to prove that the final formulas 
( [7. 71 ), ( [7. 81 ) remain true, if just to change the normalization making them non-singular for 
i G Z+: 

r(2A)r(2z)r(2z + 2A) 



Qf7v=i)(±A) 



T{2X-2iri)T{2z-2iri)T{2z + 2X + 2iri) 



X 



X 



z + X 
V 



X 
V 



9 p _P- p-'^vdz 

V 



For i G the series is terminating. In this case the above formula can be rewritten in a 
more explicit form: 



{2kri-e-n+X)dz 



19) 



k=0 



where 



9j2z + 2(£-j»0i(2A + 2(£-j» 

j=0 



9i{2z + 2X + 2{k-j)ri) 



k-l 



^2z-2{i-j)r])e,i2X-2i£-j)r]) 



X TT 

f}^ e,i2z + 2X + 2{k+j-£)7]) 



(8.10) 
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U k = oi k = i, then the second (respectively, the first) product is absent. 



The XXZ-case 

The construction of Sees. 5-7 admits several different trigonometric degenerations. Here 
we discuss the one corresponding to the XXZ spin chain with the standard 6-vertex 
i?-matrix. The L-operator is (cf. 



fS.lll 



(g27rm+jo _ g-27rm-jo j_ \ 

g27rm-jo _ g-27riM+jo J 



where u is the spectral parameter and jo, j± are operators acting on functions of a variable 
X as follows: 

jo = id, , j± = e±2--(g^e^^^^ - g-^e±^^^) . (8.12) 
Here q = e'^'^^^, and these operators obey the commutation relations of the quantum algebra 

[jo, j±] = ±27ri7j± , [j+, j_] =2(g-g ^)sinhjo. (8.13) 

The quantum transfer matrix T{u) is defined by the same formula (|2.3|) . 

The limiting transition from the elliptic formulas to the trigonometric ones consists of 
two steps: 

i) The modular transformation r — 1/t with the change of variables Zj = txj + |, 

A = TU, T] = T'j; 

ii) Tending r to 0. 

Then the L-operator defined as L^-^^^^u) = lim^^io L{X) with the understanding 

that T in the L(A) ( pj.lj ) is replaced by — 1/t, and u = A/r is kept fixed coincides with 
( p. lip . The same limiting transition in ( |7.7| ) yields the trigonometric Q-operator. 

Skipping the details (for the modular transformation of the 6'-functions and the elliptic 
gamma-function see eqs. ( |A4| ) , ( |A11| ), respectively), we present the final result. It is written 
in terms of the following special functions. The basic hypergeometric function 20i |^ is 
defined as 

^ [ J _ V MMMfc^A: .o.4^ 

201 ,q,x -2^7 ^7 , l«--L4j 

where (a; q)k = (1 — a)(l — ag)(l — ag^) ... (1 — aq^~^). Up to simple redefinitions of 
parameters, this series gives the little g-Jacobi function (when the series is terminating, it 
is the little g-Jacobi polynomial). The g-gamma- function is 
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where, for |g| < 1, (a; g)oo = liio^k^ooid] q)k- The Q-operator is given by 

ig2 [- 



Q+(n) 



-,-2mNtu 



N 

X 111201 



AT 

n 

j=i 
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p—Amu „—2£p2iTi{xj^i—Xj) 
^2g27ri(xj+i^j— 2ii) 



X 



(8.16) 



X exp{(M + £7)E^=i5x., 



As in ( [7. 81) , the formula for Q (u) differs from (|8.16|) by the changes 



u 



—U, Xi > 



and the additional cyclic shift C. The both operators Q^(m) obey the T-Q relation: 

T(m)Q(m) = sin^ 27i{u - i-f)Qiu + 7) + sin^ 27r{u + £-f)Qiu - 7) . (8.17) 

It is worth noting that the limit from (O) gives an additional operator factor which is 
omitted in (|8.16|) for it belongs to the center of the algebra generated by the shifts e^'''^""-' 
and the operator J2f=i ■ This factor has the form 

^ Tg2{-2u + 2i-f + l)Tq2{xi+i-Xi + 2i-f + l) 



n 

i=l 



rg2{xi+i-Xi-2u + l) 



where q = exp (~f^) is the modular partner of the q = exp(27ri7) (cf. [|22|). It seems so 
that the non-degenerate Sklyanin algebra knows about both halves of the modular double. 



9 Concluding remarks 

In this work, using Baxter's method of vacuum vectors, we have constructed a family of 
commuting difference operators in many variables. This family provides the Q-operator 
for the quantum integrable spin chain of XYZ-tjpe in the functional realization of the 
algebra of observables. 

The form of the obtained Q-operators suggests a simple derivation of local integrals of 
motion for the XYZ chain in infinite-dimensional representations. For instance, the nat- 
ural candidates for local Hamiltonians are (Q^(=F£?7))~^(i/(iA ^ ^ Q^(A). Their exphcit 

form is, however, not very illuminating. The hermiticity properties of these operators are 
to be further investigated. We leave this for a future work. 

Another important point remaining out of the scope of this paper is the limit to the 
XXX-case (7 ^ in ( ^.11| )). In this limit one should be able to link our results with 



the ones of the papers p, 0, [TT[, where the Q-operator for the XXX spin chain and allied 
models with rational i?-matrix was realized as an integral operator. However, it is hard to 
compare them directly since the limit is singular and, therefore, very delicate. Moreover, 
as it should be clear from the above, the rational limit is certainly not unique. Different 
rational limits of the trigonometric L-operator lead to L-operators for the XXX spin 
chains, where the sZ(2)-algebra is realized in different ways. A detailed study of rational 
limits seems to be interesting since in this case, closely related to the representation theory 
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of SL{2, M) or SU{1, 1), the result must admit a meaningful translation into the language 
of representation theory and harmonic analysis. This may give a hint of how an elliptic 
extension of the latter might look like. 

The construction of the pre-Q-operators, Ql(A) and Qk(A), goes through also for 
inhomogeneous chains. The problem, however, is the normalization: in general there is 
no point A^, such that Ql(A^,) or Qr(A*) could be explicitly inverted. Generally speaking, 
the kernel of the Q-operator can not be represented as a product of local factors, so local 
integrals of motion for inhomogeneous chains probably do not exist. In other respects, 
the pair of pre-Q-operators is a good substitute for the Q-operator. Their kernels are 
factorized though the operators do not commute. Presumably, they still obey the exchange 



commutation relations (|6.1|) . Their proof given in Sec. 6 should basically go through in the 
inhomogeneous case, too, since the key element of the proof is local in nature. 

At last, we note that the technique developed in the present paper might help to 
construct the fundamental i?-matrix in infinite-dimensional representations of the Sklyanin 
algebra. 
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Appendix A 

Theta-functions 

We use the following definition of the Jacobi ^-functions: 

Hz\r) = E exp fvrir(A; + + 2niz{k + ^] 



(Al) 

63{z\t) = ^ exp {irirk'^ + 2'Kizk^ , 



d/^{z\T) = ^ exp (ixirk'^ + 2'ni{z H — )k \ . 

fcez ^ ^ ^ 



They also can be represented as infinite products. The infinite product representation for 
the 9i{z\t) reads: 

oo 

e,{z\r) = I exp (— - n:z) (l - e'"''") (l - e2-«'=-^)^+^)) (l - e^-^'^^-^)) (A2) 
^ k=i 
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Throughout the paper we write 6a{x\T) = 6a{x). The transformation properties for shifts 
by the periods are: 

9a{x ± 1) = , 9a{x ± r) = {-lY^'^+'-'^e-^'^^^^'^9a{x) . (A3) 

Under the modular transformation r — 1/t the ^-functions behave as follows: 







r) 




-1/t) 




> 


r) 


= ,Ji/^e-^''^/^6i{z/T\ - 




03 


> 


r) 


= yV^e-^*^'/^e3(2/r| - 


1/r), 


04 


[z 


r) 




1/r). 



By 0„ we denote the space of ^-functions of order n, i.e., entire functions F{x), x G C, 
such that 

F{x + 1) = F{x) , F{x + r) = (_i)«e-™"-2™"'F(x) . (A5) 

It is easy to see that dimB„ = n. Let F{x) G 0„, then F{x) has a multiplicative 
representation of the form F{x) = cHiLi^il^ ~" ^^i), Y17=iXi = 0, where c is a constant. 
Imposing, in addition to ([ASD , the condition F{—x) = F{x), we define the space 0+ C 9„ 
of even ^-functions of order n, which plays the important role in representations of the 
Sklyanin algebra. If n is an even number, then dimB^ = + 1. 



Elliptic gamma-function 

Here we collect the main formulas on the elliptic gamma- function P, We use the 
(slightly modified) notation of |jT9[. The elliptic gamma-function is defined by the double- 
infinite product 

00 ^ _ ^2TTi{{k + l)T+{k' + l)T'-z) 

T{z\t,t')= n 1 _ ^2.^(kr+k'r'+z) (^6) 

k,k'=0 ^ ^ 

A sufficient condition for the product to be convergent is Imr > 0, Imr' > 0. We need 
the following properties of the elliptic gamma-function: 

r(2 + l|r,r') =r(^|r,r'), (AT) 

Y{z + r|r, r') = -ze-^r7^i(r')e"*"^i(2|r')r(^|r, r') , (A8) 
r(z + r'|r, r') = -ze-^r/^i(r)e"*"ei(z|r)r(z|r, r') , (A9) 



where 



oo 



r/D(r) =ei2 



— e 



2-niki 



k=l 



is the Dedekind function. Another useful property is 



T{z\T,T')T{r' -z\T,r') = . , ^"^^ / . (AlO) 
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Under the modular transformation r —^/t the eUiptic gamma-function behaves as 

(All) 



follows IS 



Y(.\r r') = e*-^(^) T{z/r \ - 1/t,t' /r) 
^''^ T{{z-T)/T'\-T/T',-l/r') 



where 



P(z) 



1 o r + r' - 1 . r2 + r'2 + 3rr' - 3r - 3r' + 1 
H : : ^ 



3rr' 2rr' 

(r + r'-l)(r + r'-rr') 
12tV 



6rr' 



(A12) 



Let us list the most frequently used formulas for r(2;) = r(2;|r, 27]). Using ( [A5| ) several 
times, we obtain: 



r(a; — 2kri) 



k-l 

l[{9i{-x + 2ri + 2jri) 

j=0 



(A13) 
(AM) 



where i? = ze^^(^+^/6)r/^(r). 

As is seen from (|A6|), the function r(z|r, 2//) has zeros at the points z = 2{k + 1)1] + 



(m + l)r + n, and simple poles at the points z = —2k'q — mr + n, where k^m run over 
non-negative integers and n over all integers. The residues of the elliptic gamma-function 
at the poles ai z = —2kri, k = 0,1,2, .. . are: 



res 



Tiz) 



where 



z=—2kri 



tq = res 



(A15) 



r(.) 



3 7ri(r+2r?)/12 



2=0 ' ' 2TTir]DiT)r]Di2ri) 



The function $^(z) ( [4. 14] ) can be extended to arbitrary complex values of i as follows: 



Elliptic hypergeometric series 

Here we follow [|1^. We define the elliptic Pochhammer symbol (the shifted elliptic facto- 
rial) by 

[x]k = [x][x + l]...[x + k-l], (A17) 
where [x] = 9\{2xrj) (cf. ( |3.5| )). By definition, the elliptic hypergeometric series is 



,...,ar+i;z\r],T) =^z^ 

fc=0 



[ai + 2k][ai]k 



r-2 



n 



[ai-am+3 + l]A: 



(A18) 
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This is an elliptic analog of the very- well-poised basic hypergeometric series [^]. The 
series is said to be balanced if z = 1 and 



r-2 



r - 5 + (r - 3)ai = 2 ^ 0^+3 



(A19) 



m=l 



For a series J2k>o Cfc of the form ( |A18|) , the balancing condition ([A19|) means that the ratio 
Cfc+i/cfc of the coefficients is an elliptic function of k. For balanced series ( |A18| ), we drop 
the argument z = 1 and the parameters i], t writing it simply as r+iW^r(«i; 0^4, • • • , ctr+i)- 
For instance, 



ioI^9(ai;a4,a5, 



aio 



k=0 



[ai-Q;m+3 + l]fc 



(A20) 



The series is called terminating if at least one of the parameters 04, ... , Or+i is equal 
to a negative integer number. In this case there is no problem of convergence. If, say 
Or+i = —n, then the series terminates at k = n. The terminating balanced hypergeometric 
series were shown [0] to possess nice modular properties. That is why they were called 
modular hypergeometric series. 

The modular hypergeometric series obey a number of spectacular identities. One of 
them is the elliptic analog of the Bailey transformation for the iqWq p7| : 



10 



Wgiai] 04, 



ioW^9(A;/?4,...,/3io) 



[ai + l]n[/3l-/?7+ !]„[/?! -/38 + l]n[ai-a7-Q8 + l]n 
[Pl + l]n[ai-a7+l]n[ai-as + l]n[Pl-P7-P8 + ^]n 



(A21) 



where aj obey the balancing condition 2 + 3ai = Z]L=i (^m+s, «io = and the param- 
eters (3j follows: 



f3i = 2ai + 1 — 04 — as — ag , 
Pi = Pi - ai + ai , 
Pi = ai, 



i = 4,5,6, 



(A22) 



10, 



This identity was proved in |T7|. 



Appendix B 

Here we present the proofs of the intertwining relations (|3.6|) and the commutation relation 

(Ol). 



For general complex values of £ the operator W£ is given by the series ( 7.13|) : 

32-(2£+l).r(2z) / Z 



We 



r(2z + 2(2£ + l)r/)* ^ ' ' ^* ^ 



The proof of ( p.6| ) is a straightforward computation. Substituting (pl^), (|El| ) into ( p.6| ) 
and comparing the coefficients of e2(^+^~^)^^^, A; > 0, we get, after cancelling out common 
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factors: 



+ 9i{2{i-k + l)r])9^i{-2z+2{i+l)r])9i{2z+4kr])9i{2z+2{i+k + l)r]) 



9i{2{i+k + l)r])9^i{2z + 2{2k-£-l)r])9i{2z)9i{2z + 2{i+k + l)7]) 



+ 9i{2{i-k+l)ri)9c^i{-2z-2{2k+£+l)ri)9i{2z)9i{2z-2{i-k+l)ri) . 

All the four terms, as functions of 2z, have the same monodromy properties. They are 
^-functions of 3-d order. The linear space of such functions is 3-dimensional, so it is 
enough to verify the equality at 3 points, say z = 0, z = —2kri, and z = {£ — k + 1)7]. 
The coefficients of the highest shift e'^^^^^^'^^^ in both sides of ( p. 61 ) are equal, too, so the 
intertwining relation is proved. 



To prove (|7.15|) , we use the representation ( p.2|) . For any two adjacent sites of the 

12 12 

chain (say 1 and 2) we define the c-number matrix M=M = V {\^, zi)V{\+, Z2) 

explicitly given by 



12 



M {X,i)- 



1 



9i {zi + Z2- 2ir])9i (^1 -Z2 + 2X) -9^ (^i - ^2 - 2ir])9i{zi + Z2 + 2X) 



^1(2^2) \ 0^^^^_^^^2ir])9i{zi+Z2-2X) -9i{zi + Z2 + 2ir])9i{zi-Z2-2X) 



Using the cyclic property of the trace, one can rewrite (|2.3| ) as follows: 



T(A) = :tr 



12 / e''^^ \ 23 / e''^^ 



M 



where di = dz^. Computing the trace, we obtain: 



ATI / QVdi Q 



.. M 



e"''^! 



12 23 



Afl 



T(A) = Mi.i^Mi^r, ■ ■ ■ Mi^i, e 

il,...,iiv=l,2 



(3-2ji)»)9ig(3-2j2)»)92 g(3-2«jv)»79jv 



fB2) 



Recall that (pi{zi, Zi+i) = ^e{zi — Zi^i)^£{zi + Zi^i) . The commutation relation (|7.15| ) 



N N 



i=l i=l 

follows from (IB 21) and the local commutation rule 



e^^-^'''^''^'^'^-^'^^''^^^Azi,Z2) 



9i{zi + {-ly'+'^z2-2{-ly'{^ + l)r]) 



■(pe[zi,Z2)e 



(3-2n)r;ai4{3-2j2)»?92 



, Zi,Z2 = 1,2 



9,{zi + {-iy^+'^Z2 + 2{-iy^ir]) 
which is easy to verify. Writing out the entries of the matrix M in the form 

M^k= 9i{z, + {-iy+''z2 + 2{-iye7])Qk , t,k = l,2, (B3) 

(the second factor stands for something that does not depend on i), and subsequently 
moving the factors (p£{zi, Zi^i) from the right to the left through each term of (p^) , we 
arrive at ( |7.15| ). 



29 



References 



[1] R.Baxter, Partition function of the eight-vertex lattice model, Ann. Phys. 70 (1972) 
193-228; R.Baxter, Eight-vertex model in lattice statistics and one- dimensional 
anisotropic Heisenberg chain. I, II, III, Ann. Phys. 76 (1973) 1-24, 25-47, 48-71 

[2] L.D.Faddeev and L.A.Takhtadzhan, The quantum method of the inverse problem and 
the Heisenberg XYZ model (English translation), Russian Math. Surveys 34:5 (1979) 
11-68 

[3] I.G.Korepanov, Vacuum curves of C-operators connected with the six-vertex model 
and construction of TZ-operators, Chelyabinsk Polytechn. Inst., Deposited in the 
VINITI 02.04.86, N 2271-B86; Algebra i Analiz 6:2 (1994) 176-194; V.V.Bazhanov 
and Yu.G.Stroganov, Chiral Potts model as a descendant of the six-vertex model, J. 
Stat. Phys. 59 (1990) 799-817 

[4] M.Batchelor, R.Baxter, M.O'Rourke and C.Yung, Exact solution and interfacial ten- 
sion of the six-vertex model with anti-periodic boundary conditions, J. Phys. A 28 
(1995) 2759-2770; M.Batchelor and C.Yung, Exact solution for the spins XX Z quan- 
tum chain with non-diagonal twists, Nucl. Phys. B446 (1995) 461-484 

[5] V.Pasquier and M.Gaudin, The periodic Toda chain and a matrix generalization of 
the Bess el function recursion relations, J. Phys. A 28 (1992) 5243-5252 



[6] S.Derkachov, Baxter's Q-operator for the homogeneous XXX spin chain, |solv- 



I int/99020T5 

[7] G.Pronko, On the Baxter's Q-operator for the XXX spin chain, Piep-th/9908179| 

[8] V.Bazhanov, S.Lukyanov and A.Zamolodchikov, Integrable structure of conformal 
field theory 11. Q-operator and DDV equation, Commun. Math. Phys. 190 (1997) 
247-278; Integrable structure of conformal field theory III. The Yang-Baxter relation, 
Commun. Math. Phys. 200 (1999) 297-324 

[9] A.Antonov and B.Feigin, Quantum group representations and Baxter equation, Phys. 
Lett. B392 (1997) 115-122 

[10] V.Kuznetsov and E.Sklyanin, On Bdcklund transformations for many-body systems, 
J. Phys. A 31 (1998) 2241-2251 

[11] V.Kuznetsov, M.Salerno and E.Sklyanin, Quantum Bdcklund transformation for the 
integrable BST model, |solv-int / 990800^ 

[12] T.Takebe, Generalized Bethe ansatz with the general spin representation of the 
Sklyanin algebra, J. Phys. A 25 (1992) 1071-1083; Bethe ansatz for higher spin eight 
vertex models, J. Phys. A 28 (1995) 6675-6706 

[13] V.Bazhanov and N.Reshetikhin, Critical RSOS models and conformal field theory. 
Int. J. Mod. Phys. A4 (1989) 115-142 

[14] I.Krichever, O.Lipan, P.Wiegmann and A.Zabrodin, Quantum integrable models and 
discrete classical Hirota equations, Commun. Math. Phys. 188 (1997) 267-304 



30 



[15] E.K.Sklyanin, On some algebraic structures related to the Yang-Baxter equation, Funk. 
Anal, i ego Pril. 16:4 (1982) 27-34 

[16] E.K.Sklyanin, On some algebraic structures related to the Yang-Baxter equation. Rep- 
resentations of the quantum algebra, Funk. Anal i ego Pril. 17:4 (1983) 34-48 

[17] I.Frenkel and V.Turaev, Elliptic solutions of the Yang-Baxter equation and mod- 
ular hypergeometric functions, in: The Arnold- Gelf and Mathematical Seminars, 
Birkhauser, Boston, MA, 1997, 171-204 

[18] S.N.M.Ruijsenaars, First order analytic difference equations and integrable quantum 
systems, J. Math. Phys. 38 (1997) 1069-1146 

[19] G.Felder and A.Varchenko, The elliptic gamma function and SL(3,Z) x Z^, 
|math.QA/9907061| 

[20] A.Zabrodin, On the spectral curve of the difference Lame operator. Int. Math. Research 
Notices No. 11 (1999) 589-614 

[21] I.Krichever, Baxter's equations and algebraic geometry. Funk. Anal, i ego Pril. 15:2 
(1981) 22-35 



[22] L.D.Faddeev, Modular double of quantum group, |math.QA/9912078| 



[23] L.N.Lipatov, High energy asymptotics of multi-colour QCD and exactly solvable lattice 
models, JETP Lett. 59 (1994) 571-574 

[24] L.Faddeev and G.Korchemsky, High energy QCD as a completely integrable model, 
Phys. Lett. B342 (1995) 311-322 

[25] I.Krichever and A.Zabrodin, Spin generalization of the Ruijsenaars-Schneider model, 
non-abelian 2D Toda chain and representations of Sklyanin algebra, Usp. Mat. Nauk, 
50:6 (1995) 3-56, |hep-th/9505039| 



[26] G.Felder and A.Varchenko, Algebraic integrability of the two-body Ruijsenaars opera- 
tor. Funk. Anal, i ego Pril. 32 (1998) 8-25, |q-alg/961002^ 

[27] G. Gasper and M.Rahman, Basic hypergeometric series. Encyclopedia of Mathematics 
and its Applications, Volume 35, Cambridge University Press, 1990 



31 



